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and report subsequent drastic boost of the effective diffusivity for weak advection. 
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Disorder in operation conditions of a dynamic system is known to be able to play not 
only trivial destructive role, distorting the system behavior, but also constructive one, 
inducing certain degree of order and leading to various non-trivial effects: Anderson 
localization [1, 2, 3, 4], stochastic [5] and coherence resonances [6], noise-induced 
synchronization [7], etc. One of the most distinguished and fair effects is the Anderson 
localization (AL), which is the localization of states in spatially extended linear systems 
subject to a frozen parametric disorder (random spacial inhomogeneity of parameters). 
AL was first discovered and discussed for quantum systems [1]. Later on, investigations 
were extended to diverse branches of semiclassical and classical physics: wave optics [3], 
acoustics [4], etc. The phenomenon was comprehensively studied and well understood 
mathematically for the Schrodinger equation and related mathematical models [2, 8, 9]. 
The role of nonlinearity in these models was addressed in literature as well (e.g., [9, 10]). 

While extensively studied for conservative media, the localization phenomenon did 
not receive comparable attention for active/dissipative ones, like in problems of thermal 
convection or reaction-diffusion. The main reason is that the physical interpretation 
of formal solutions to the Schrodinger equation is essentially different from that of 
governing equations for active/dissipative media; therefore, the theory of AL may be 
extended to the latter only under certain strong restrictions (see [11, 12] for reference). 
Nevertheless, effects similar to AL can be observed in fluid dynamical systems [11, 12]. 
In [12] we addressed the problem where localized thermoconvective currents excited 
in a horizontal porous layer under frozen parametric disorder (spacial inhomogeneity 
of the macroscopic permeability, the heat diffusivity, etc.) drastically influence the 
process of transport of a passive scalar (e.g., a pollutant) along the layer. Below the 
threshold of instability of the disorder-free system, the effective diffusivity quantifying 
this transport has been found to be faithfully determined by localization properties of 
patterns. Meanwhile, in [11] these properties have been revealed to be greatly affected 
by a weak imposed longitudinal advection. Hence, one can expect a weak advection 
to lead to a significant enhancement of the effective diffusivity of a nearly indiffusive 
pollutant. Treatment of this effect is the subject of the present paper. 

The paper is organized as follows. In section 1 we formulate the specific 
physical problem that we deal with and introduce the relevant mathematical model; 
in particular, we discuss disorder-induced excitation of localized currents below the 
instability threshold of the disorder-free system and advectional derealization of these 
patterns. Section 2 presents the results of a numerical simulation and calculation of 
the effective diffusivity. In section 3 we develop an analytical theory for the effective 
diffusivity in the presence of an imposed advection. Section 4 ends the paper with 
conclusions. 

1. Problem formulation and current state of research 

In nature and technology, a broad variety of active media where pattern selection occurs 
is governed by Kuramoto-Sivashinsky type equations. In the presence of an imposed 
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advectional transport u in the x-direction the modified Kuramoto-Sivashinsky equation 
reads 



This equation describes two-dimensional large-scale natural thermal convection in a 
horizontal fluid layer heated from below [13, 14] and is still valid for a turbulent fluid [15], 
a binary mixture at small Lewis number [16], a porous layer saturated with a fluid [17], 
etc. In these fluid dynamical systems, except for the turbulent one [15], the plates 
bounding the layer should be nearly thermally insulating (in comparison to the fluid) 
for a large-scale convection to arise. In the problems mentioned, equation (1) governs 
evolution of temperature perturbations 9 which are nearly uniform along the vertical 
coordinate z and determine fluid currents. 

The origin of such a frequent occurrence of equation (1) is its general validity, which 
may be argued as follows. Basic laws in physics are conservation ones. This often results 
in final governing equations having the form d t [quantity] + V- [flux of quantity] = 0. Such 
conservation laws lead to Kuramoto-Sivashinsky type equations. While the original 
Kuramoto-Sivashinsky equation has a quadratic nonlinear term (cf [18]), this term 
should be replaced by a cubic one for the systems with the sign inversion symmetry 
of the fields, which is widespread in nature, or for description of a spatiotemporal 
modulation of an oscillatory mode. Thus the governing equation takes the form (1). 
On these grounds, we state that equation (1) describes pattern formation in a broad 
variety of physical systems. 

We restrict this paper to the case of convection in a porous medium [17, 12] 
for the sake of definiteness; nonetheless, most of our results may be extended in a 
straightforward manner to the other physical systems mentioned. Equation (1) is already 
dimensionless and below we introduce all parameters and variables in appropriate 
dimensionless forms. 

In the large-scale (or long-wavelength) approximation, which we use, the 
characteristic horizontal scales are assumed to be large against the layer height h. 
In equation (1), q(x) represents the local supercriticality: (21/2)h 2 q(x) is the sum of 
relative deviations of the heating intensity and of the macroscopic properties of the 
porous matrix (porosity, permeability, heat diffusivity, etc.) from the critical values for 
the spatially homogeneous case [17]. For positive spatially uniform q, convection sets up, 
while for negative q, all the temperature perturbations decay. In a porous medium [17], 
the macroscopic fluid velocity field is 



where ip(x,t) = 6 x (x,t) is the stream function amplitude, and the reference frame 
is such that z = and z = h are the lower and upper boundaries of the layer, 
respectively [figure 1(b)]. Though the temperature perturbations obey equation (1) 
for diverse convective systems, function f(z), which determines the relation between 
the flow pattern and the temperature perturbation, is specific to each case. Notice that 



0(x,t) 



(u 9{x, t) + 9 xxx {x, t) + q(x) 9 x (x, t) - (9 x (x, t)f) 



(1) 



oz ox 
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Figure 1. (a) Establishing steady solutions to equation (1) are sets of exponentially 
localized patterns [shown for one and the same realization of random inhomogeneity 
£(x) and e = 1; q(x) is represented by q 7I (x) = n^ 1 l{ x ')dx']- 
(b) The stream lines corresponding to the solutions in graph (a) are plotted for the 
case of convection in a porous layer [see equation (2) for an exact relation] . 

u is not presented in expression (2) owing to its smallness in comparison to the excited 
convective currents v. The impact of a weak imposed advective flow on the evolution 
of temperature perturbations is caused by its symmetry properties: the gross advective 
flux through the vertical cross-section is u, while the convective flow v possesses zero 
gross flux and, therefore, yields a less effective heat transfer along the layer [17]. 
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Figure 2. Spectrum of localization exponents of time-independent patterns for e = 1, 
go = — 1 derived from equation (4) [11]. 

Although equation (1) is valid for a large-scale inhomogeneity q(x), which means 
ft-kzl/kl ^ 1) one can set such a hierarchy of small parameters, namely h <C 
(h\q x \/\q\) 2 <C 1, that a frozen random inhomogeneity may be represented by white 
Gaussian noise (,(x): 

q(x) =q + (£(*)) = 0, (tfsKOO) = ^(x ~ A 

where e 2 is the disorder intensity and go is the mean supercriticality (i.e. departure from 
the instability threshold of the disorder-free system). Numerical simulation reveals only 
time-independent solutions to establish in (1) with u — and such q(x) [11]; for a small 
non-zero u, stable oscillatory regimes are of low probability by continuity. 
In the stationary case for u — the linearized form of equation (1), i.e., 

-Oxxx(x) - £(x) 9 x (x) = q 9 x (x) , (3) 

is a stationary Schrodinger equation for ip = 9 X with q instead of the state energy 
and —£(x) instead of the potential. Therefore, similarly to the case of the Schrodinger 
equation (see [2, 8, 9]), all the solutions ip(x) to the stationary linearized equation (1) 
are spatially localized for arbitrary q ; asymptotically, 

ip(x) oc exp(— 7|x|), 

where 7 is the localization exponent. Such a localization can be readily seen for the 
solution to the nonlinear problem (1) in figure 1(a) for go — —2.5, u — 0, which is a 
solitary vortex. 

For u 7^ the solitary patterns are still exponentially localized [figure 1(a), u = 0.3]. 
However, their localization properties change drastically because instead of the second- 
order linear ODE with respect to ip — 9 X , equation (3), one finds a third-order equation: 

-u9(x) - 9 xxx (x) - £(x) 9 x (x) = q 9 x (x) . (4) 

The two symmetric modes 9 oc exp(±7a;) and trivial solution 9 = const (ip = 0) turn 
into three modes 9 oc exp(7jx), 73 < < 72 < 71, of equation (4) with u > 0, q < 
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(see sample spectrum of 7$ in figure 2; cf [11] for details). Specifically, 71-mode is the 
successor of +7, 72-mode is the one of the trivial homogeneous mode, and 73-mode is 
that of —7. Thus, the upstream flank of the localized pattern is now composed by two 
modes decaying in the distance from the pattern: 

0(x) « 0i(a;)e 7ia; + Q 2 {x)e~ /2X or ij>(x) « 7i9i(x)e 7ia; + 7 2 6 2 (a;)e 72X , 

where functions Qi(x) and 6«(x) neither grow nor decay over large distances. The 71- 
mode, which disappears for u — 0, i.e. 71 = 0, decays slowly for a small finite u, prevails 
over the 72-mode decaying rapidly, and, thus, determines the upstream localization 
properties of the pattern. The upstream localization length 1/7! can become remarkably 
large leading to upstream derealization of patterns, which can be seen in figure 1. 

One should keep in mind, that consideration of solitary patterns makes sense where 
such patterns can be distinguished, i.e., are sparse enough in space. This is the case of 
negative qo < —1- In figure 1, for a sample realization of £(x), one can see that localized 
patterns can be discriminated for e — 1, q — —1.5 and the localization properties are 
very well pronounced for q = —2.5. 

Here we would like to emphasize the fact of existence of convective currents below 
the instability threshold of the disorder-free system. These currents considerably and 
nontrivially affect transport of a pollutant (or other passive scalar), especially when its 
molecular diffusivity is small in comparison to the thermal one, which is quite typical in 
nature (for instance, at standard conditions the molecular diffusivity of NaCl in water is 
1.1 ■ 10 _9 m 2 /s against the heat diffusivity of water which is 1.3 • 10 _7 m 2 /s). Transport 
of a nearly indiffusive passive scalar, quantified by the effective (or eddy) diffusivity 
coefficient, is the object of our research, as a "substance" which is essentially influenced 
by these localized currents and, thus, provides an opportunity to observe manifestation 
of disorder-induced phenomena discussed in [11]. 

In [12] we studied the problem for the case of no advection (u = 0) and calculated 
(both numerically and analytically) the enhancement of the effective diffusivity by 
disorder-induced currents; this enhancement is especially strong for low molecular 
diffusivity deep below the instability threshold of the disorder-free system (see figure 3). 
In this paper we address the role of an imposed advection in this problem. The interest 
to advection is provoked by its dramatic influence on localization properties, i.e., the 
upstream derealization of convective currents that is described above. We expect 
this derealization to result in a giant increase of the effective diffusivity for a nearly 
indiffusive pollutant and, particularly, in the lowering of the mean supercriticality (go) 
value at which the transition from sets of localized convective currents to an almost 
everywhere intense 'global' flow occurs. 

2. Effective diffusivity 

In this section we describe the transport of a passive pollutant by a steady convective 
flow (2); "passive" means that the flow is not influenced by the pollutant. The 
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assumption of passiveness is practically relevant because (biologically/chemichally) 
significant concentrations of a pollutant can be very small and mechanically negligible. 
The flux j of the pollutant concentration C is 

J=vC-DVC, (5) 

where the first term describes the convective transport, the second one represents 
the molecular diffusion, and D is the molecular diffusivity. The establishing time- 
independent distributions of the pollutant obey 

V-J=0. (6) 

Equation (6) [with account for (2)] yields a distribution of C which is uniform along z 
and is determined by 

dC(x) J 



dx ~ „ 21^ 2 (x) ' 
D + 



(7) 



2h 2 D 

where J is the constant pollutant flux along the layer. Detailed derivation of equation 
(7) for u = can be found in [12] where it was performed in the spirit of the standard 
multiscale method (interested readers can consult, e.g., [19, 20]). Remarkably, advection 
velocity u is not presented in the last equation: its direct contribution to convective 
currents transferring the pollutant is small in comparison to the one of excited convective 
currents. Instead, it influences the heat transfer and, consequently, excited flows, 
drastically changing properties of the field ip(x). Notice that, for the other convective 
systems which we mentioned in section 1, the result differs only in the factor ahead of 
^ 2 /D. 

Thus we come to introducing the effective diffusivity for the system under 
consideration (general ideas on the effective diffusivity in systems with irregular currents 
can be found, e.g., in [21, 20]). Let us consider the domain x G [0,L] with the imposed 
concentration difference SC at the ends. Then the establishing pollutant flux J is defined 
by the integral [cf (7)] 

. r dx 
J 21^ 2 (x) ' 

For a lengthy domain the specific realization of becomes insignificant: 



SC = -J 




Hence, 



SC 

J=-a-, 

which means that a can be treated as an effective diffusivity. 
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Figure 3. Dependencies of effective diffusivity a on mean supercriticality 5o in the 
presence of an imposed advection (u = 0.2) and without it (u = 0). The bold black 
solid line represents the analytical dependence for u = 0, and the dashed line represents 
the one for u — 0.2 (see section 3). 



The effective diffusivity 

~(("^)T 

turns into D for vanishing convective flow. For small D the regions of the layer where 
the flow is damped, ip <C 1, make large contribution to the mean value appearing in (8) 
and diminish a, thus, leading to the locking of the spreading of the pollutant. 

The disorder strength e 2 can be excluded from equations by the appropriate 
rescaling of parameters and fields. As a consequence, the results on the effective 
diffusivity can be comprehensively presented in terms of D, a, %, and u: 

D=J^e^D, a = ^e^ 3 ha, q = , u = - . 
V 21 ' V 21 yu e 4 / 3 ' e 2 

Figure 3 provides calculated dependencies of effective diffusivity a on q for moderate 

and small values of molecular diffusivity D. Concerning these dependencies the following 

is worth noticing: 

(a) For small D a quite sharp transition of effective diffusivity a between moderate 
values and ones comparable with D occurs near q = (note logarithmic scale of the 
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vertical axis), suggesting the transition from an almost everywhere intense 'global' flow 
to a set of localized currents to take place. 

(b) In the presence of a weak imposed advection, u = 0.2, the transition to 'global' flow 
occurs at the value of % which is considerably lower than that without advection. 

(c) Below the instability threshold of the disorder-free system, where only sparse 
localized currents are excited, the effective diffusion can be significantly enhanced by 
these currents. 

(d) The disorder-induced enhancement of the effective diffusivity is especially drastic 
in the presence of an imposed advection; e.g., for D = 10" 4 , q = —1, the effective 
diffusivity is increased by one order of magnitude compared to the molecular diffusivity 
without advection (u = 0) and by two orders of magnitude for u = 0.2. 



3. Analytical theory 

3.1. Transport through time-independent patterns 

The effective diffusivity can be analytically evaluated for a small molecular diffusivity 
(D <C 1) and sparse domains of excitation of convective currents (the spacial density of 
the excitation domains v <C 1). In [12] it was evaluated for the case of no advection, 

az=o~D^iy, (9) 
where one can use the asymptotic expressions for the density of the excitation domains 

1 ( 1-95 qj \ nn . 



4 v / r95Te 2 / 3 |g | 



and 



2 / , 1 



= 0) = : ? [ |(/o|" - -|g | + ■ ■ ■ 



which are valid for % < — 1. The latter expression is known from the classical theory 
of AL (e.g., see [8, 9]). In the following we advance the evaluation procedure realized 
in [12] in order to account for the asymmetry between up- and downstream localization 
exponents. 

Now we calculate the average 

Due to ergodicity, this average over x for a given realization of £(x) coincides with the 
average over realizations of £(x) at a certain point Xq. We set the origin of the x-axis 
at x Q and find o~ x = ((D + ^ 2 (0)/5) -1 )^. 

When the two nearest to the origin excitation domains are distant and localized 
near x\ > and x 2 < (see figure 4), 

^(0) w 7i^ 7lXl + l20 2 e~ l2Xl + 7 3# 3 e- 73 " 2 , (11) 
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Figure 4. Sketch of two localized flow patterns being nearest to the origin. 



where 0i j2 and 9 3 characterize the amplitude of temperature perturbation modes excited 
around X\ and x 2 , respectively. For small D and density u, the contribution of the 
excitation domains to a -1 is negligible against that of the extensive regions where flow 
is weak, but is still larger than D. Therefore, one may be not very subtle with 
"cores" of excitation domains and may utilize expression (11) even for small x\ i2 '- 



D + ?P 2 (0)/D 



dxi / dx 



p(x 1 )p(x 2 ) 




D + D-^Ae-^ 1 + <y 2 9 2 e-^ + 7 3 3 e'w a:2 ) 2 

where p(x\) [p{xo)] is the density of the probability to observe the nearest right [left] 
excitation domain at +xi [—x 2 ]. For probability distribution P(xi > x), one finds 
P(xi > x + dx) = P(xi > x) (1 — vdx), i.e., (d/dx)P(xi > x) = —vP{x\ > x). Hence, 
P(xi>x) = e~ ux , and probability density p(x) = \(d/dx)P(xi>x)\ = ue~ ux . As regards 
averaging over it is important that the multiplication of 9{ by factor F is effectively 
equivalent to the shift of the excitation domain by |7i| -1 InF, which is insignificant for 
F ~ 1 in the limit case that we consider. Hence, one can assume 9i = ±1 (the topological 
difference between different combinations of signs of 9{ is not to be neglected) and rewrite 
equation (12) as 



1 

a 



dx\ J dx 2 v e 
ii 



2 -u(x 1 +x 2 ) 



1 



+ 



D + D- 1 (7ie~T lXl + ^ 2 e~ l2Xl + ^e lzX2 ) 2 
1 



D + D~ 1 (7ie~i' ia:i + ^ 2 e~ l2Xl - ^e 1 ^ 2 ) 2 

For u/ji 1 and 72-mode dominating over 7!-mode (that is the case in figure 1), the 
last formula yields 



a 



D 



72 



173 



(13) 



1731 

DJ V D 

Here we assume that advection is weak and suppresses thermal convection in a negligible 
fraction of the excitation centers and the asymptotic expression (10) is still valid. 

Noticeable difference between equation (13) for u = 0, i.e., 72 = —73 = 7, 
and equation (9) is actually insignificant up to our approximations, because moderate 
number 7/2 risen to small power 2z//7, (7 /2) 2u ^ which is the ratio of these equations, 
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Figure 5. The theoretical dependence of effective diffusivity a on mean supercriticality 
§o and molecular diffusivity D in the presence of an imposed advection, u = 0.2 (black 
"wireframe" surface), and without it, u — (gray surface). 

approximately equals l.f For instance, in figure 3, the curves given by analytic 
expressions (9) and (13) with 72 = —73 = 7 are visually undistinguishable. 

For small finite u, smallness of 72 in equation (13) gives rise to a significant 
enhancement of effective diffusivity a, which is in agreement with the results of numerical 
simulation presented in figure 3. Figure 5 shows that for u = 0.2 the effective diffusivity 
in the presence of advection is always stronger than without it; the larger the difference 
between the effective and the molecular diffusivity the stronger advectional enhancement 
of the effective diffusivity is. For instance, for D m 10~ 4 , qo — — 1 the effective diffusivity 
in the presence of advection u = 0.2 is by factor 10 larger than without advection, and 
this factor grows as D decreases. 

Noteworthy, for u = 0.2 expression (13) provides slightly overestimated value of the 
effective diffusivity while for u = the analytical estimation is accurate. The inaccuracy 
appears because in our analytical theory we ignore three factors: (a) decrease of the 
spatial density of the excitation centers owing to advectional suppression (washing- 
out) of weak excitation centers; (b) for small u the rapidly decaying upstream 71- 
mode is significant because of the smallness of the slowly decaying 72-mode; and 
(c) as the advection strengthens, currents in some excitation domains disappear via 
a Hopf bifurcation [11] and thus there is non-zero probability to observe oscillatory 
flows for small u even though there is no stable time- dependent solutions for u — 0. 
Unfortunately, inaccuracies caused by these three assumptions can not be minimized 

I Equation (9) is more accurate than equation (13), because the case of u = is simpler than that we 
consider here and admits analytical evaluation of integrals with a fewer number of approximations. 




Figure 6. (a) Sample of an oscillatory flow appearing for u = 0.2, e = 1: at 
qo = —1.5 the pattern is time-independent, further, as §o decreases, it turns oscillatory 
(§o = —1.6, shown at a certain moment of time) and then decays for go ^ —1.7. 
(b) The effective diffusivity coefficients calculated over a short domain (length L = 
100), with sample t;(x) and the patterns plotted in graph (a), are compared with the 
one in the limit of an infinite domain (L — 10 5 ). 



simultaneously: the first and third assumptions require u — > 0, while the second one 
needs u to be small but finite. 

3.2. Discussion of transport through oscillatory patterns 

Oscillatory localized patterns discovered in the dynamic system (1) for non-zero u 
(figure 6(a); see [11] for details) are statistically improbable and rare when u is small. 
Notably, their relative contribution to the effective diffusivity is much larger than 
their fraction among the excited localized patterns. In figure 6(b) one can see the 
soaring of the effective diffusivity along a finite region as a localized pattern turns 
oscillatory (g = —1.6) before disappearing (q < —1.7). Figure 6(a) reveals the origin 
of this soaring: the oscillatory pattern is not so well localized as the time-independent 
one. Indeed, the localization properties of the oscillatory pattern of frequency to are 
determined by the following linearization of equation (1): 

iu6 = -(u6 + 9 XXX + (g + Z{x))B x ) x . (14) 

In contrast to (4), this is already a 4th-order differential equation, which yields four 
localization exponents. The newly appeared 4th mode possesses 7 oc u oc u, i.e., decays 
slowly, and contributes to the downstream flank of localized patterns (evidence of these 
facts is beyond the scope of this paper and will be presented elsewhere). As well as the 
72-mode results in upstream derealization of time-independent patterns, the new mode 
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leads to downstream derealization of oscillatory patterns, which appear to be weakly 
localized both up- and downstream, as one sees this in figure 6(a). 

Nevertheless, owing to the smallness of the fraction of the oscillatory patterns 
among all the localized patterns at small u, their contribution to the effective diffusivity 
over large domains is still negligible. This is additionally confirmed by the accuracy of 
our analytical theory disregarding oscillatory currents [equation (13)]. Meanwhile, the 
analytical theory accounting for the oscillatory patterns should involve the distribution 
of frequencies of excited patterns, which are to be determined only from the nonlinear 
problem (1): this is not an analytically solvable problem. 

4. Conclusion 

We have studied the transport of a pollutant in a horizontal fluid layer by spatially 
localized two-dimensional thermoconvective currents appearing under frozen parametric 
disorder in the presence of an imposed longitudinal advection. Though we have 
considered the specific physical system, a horizontal porous layer saturated with a fluid 
and confined between two nearly thermally insulating plates, our results can be in a 
straightforward manner extended to a broad variety of fluid dynamical systems (like 
ones studied in [13, 14, 15, 16]). We have calculated numerically the dependence of the 
effective diffusivity on the molecular one and the mean supercriticality for a non-zero 
advection strength (see figure 3). The results reveal that advectional derealization of 
convective currents greatly assists transfer of a nearly indiffusive pollutant (D < 1) 
below the instability threshold of the disorder-free system: the effective diffusivity can 
become by several orders of magnitude larger in comparison to that without advection. 

The analytical theory focusing on advectional derealization of localized current 
patterns yields results which are in a fair agreement with the results of numerical 
simulation. This correspondence confirms our treatment of importance of disorder- 
induced patterns and their localization properties in active/dissipavite media, which 
provoked works [11, 12]. 
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